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Problem  Statement 


The  objective  of  this  work  was  to  develop  and  demonstrate  a  combined 
Monte  Carlo/finite-difference  procedure  for  calculating  the  properties  of 
two-dimensional  combustor  flows. 

It  is  a  formidable  theoretical  problem  to  calculate  the  turbulent  reacting 
flow  properties  in  a  gas  turbine  combustor.  Added  to  the  well-known  problems 
of  calculating  any  turbulent  flow,  are  the  difficulties  caused  by  the  highly 
non-linear  combustion  reactions.  An  approach  that  goes  a  long  way  to  over¬ 
coming  the  problems  caused  by  the  non-linear  reactions  is  to  describe  the 
turbulent  reacting  flow  field  in  terms  of  the  joint  probability  density 
functions  (pdf)  of  the  species  mass  fractions.  The  transport  equation  for 
the  joint  pdf  of  the  species  contains  an  exact  term  pertaining  to  reaction. 

That  is,  irrespective  of  the  complexity  and  non-linearity  of  the  reaction 
scheme,  in  the  joint  pdf  equation  the  reaction  term  does  not  have  to  be 
modelled.  For  this  reason,  several  authors  have  used  this  equation  to  calculate 
turbulent  reactive  flows:  Dopazo  (1975),  Janicka,  Kolbe  and  Kollmann  (1979) 
and  Pope  (1981a) . 

By  standard  techniques  (e.g.  finite-differences),  the  joint  pdf  equation 
would  be  extremely  expensive  to  solve  numerically  (if,  ineeed,  it  were 
possible),  except  in  simple  cases.  However,  by  using  a  Monte  Carlo  method, 
the  equation  can  be  solved  at  a  reasonable  cost  for  situations  of  practical 
complexity.  The  basic  Monte  Carlo  method  is  described  by  Pope  (1981b). 

The  Monte  Carlo  solution  of  the  joint  pdf  equation  deals  only  with  the 
species  fields.  The  mean  velocity  and  turbulence  fields  are  calculated  by  a 
standard  finite  difference  procedure  using  the  k-c  turbulence  model. 
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Results 

In  the  Monte  Carlo  method  described  by  Pope  (1981b),  the  joint  pdf  is 
represented  by  ensembles  of  elements  at  discrete  grid  nodes.  There  is  both 
a  truncation  error  and  a  statistical  error  associated  with  the  discretization 
of  physical  space.  An  improved  method  has  been  developed  in  which  the 
elements  are  distributed  throughout  the  flow  field  —  not  just  at  discrete 
nodes.  This  method  resembles  a  Lagrangian  particle  tracking  procedure  with 
a  random  walk  to  simulate  turbulent  diffusion.  The  numerical  error  is 
greatly  reduced,  and  hence  the  computational  economy  is  improved.  The  method 
is  described  by  Pope  (1982a) . 

This  method  has  been  inplemented  for  plane,  constant-density  flows.  The 
two  major  tasks  were  to  treat  the  boundaries  correctly  and  to  develop  the 
necessary  least-squares  cubic  spline  procedures.  The  cubic  splines  are  used 
to  extract  information  about  mean  quantities  from  the  elements  randomly 
distributed  through  the  flow  field.  At  the  end  of  the  contract,  the  code 
had  been  developed  and  testing  was  under  way. 

One  of  the  components  of  the  Monte  Carlo  method  is  a  model  of  turbulent 
mixing.  Most  investigators  use  Curl's  (1963)  model  or  a  variant  of  it 
(Dopazo,  1979;  Janicka,  Kolbe  and  Kollmann,  1979).  In  the  course  of  this 
investigation  it  was  found  that  these  models  have  serious  defects.  Specifi¬ 
cally,  when  the  models  are  applied  to  the  simple  case  of  decaying  fluctuations 
they  should  predict  a  Gaussian  pdf.  However,  the  predicted  pdf  is  far  from 
Gaussian  and,  in  fact,  its  fourth  and  higher  even  moments  are  infinite. 

Because  of  this  observed  defect,  a  new  class  of  mixing  models  has  been  developed 
that  overcomes  this  problem  (Pope,  1982b).  The  improved  models  result  in 
approximately  Guassian  pdf's,  and  a  consequent  gain  in  accuracy  is  expected 
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"An  Improved  Turbulent  Mixing  Model,"  by  S.  B.  Pope,  to  be  published  in 
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"PDF  Methods  in  Turbulent  Reactive  Flows,"  by  S. 
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